Computational Design of Chemical Nanosensors: Metal Doped Carbon Nanotubes 
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We use computational screening to systematically investigate the use of transition metal doped carbon nan- 
otubes for chemical gas sensing. For a set of relevant target molecules (CO, NH3, H2S) and the main components 
of air (N2, O2, HtO), we calculate the binding energy and change in conductance upon adsorption on a metal 
atom occupying a vacancy of a (6,6) carbon nanotube. Based on these descriptors, we identify the most promis- 
ing dopant candidates for detection of a given target molecule. From the fractional coverage of the metal sites 
in thermal equilibrium with air, we estimate the change in the nanotube resistance per doping site as a function 
of the target molecule concentration assuming charge transport in the diffusive regime. Our analysis points to 
Ni-doped nanotubes as candidates for CO sensors working under typical atmospheric conditions. 



PACS numbers: 73.63.-b, 68.43.-h, 73.50.Lw 

The ability to detect small concentrations of specific chem- 
ical species is fundamental for a variety of industrial and sci- 
entific processes as well as for medical applications and en- 
vironmental monitoring [1]. In general, nanostructured mate- 
rials should be well suited for sensor applications because of 
their large surface to volume ratio which makes them sensi- 
tive to molecular adsorption. Specifically, carbon nanotubes 
(CNT) [2] have been shown to work remarkably well as de- 
tectors of small gas molecules. This has been demonstrated 
both for individual CNTs [3-8] as well as for CNT networks 
[9, 10]. 

Pristine CNTs are known to be chemically inert - a prop- 
erty closely related to their high stability. As a consequence, 
only radicals bind strong enough to the CNT to notably affect 
its electrical properties [2, 5, 11-13]. To make CNTs attrac- 
tive for sensor applications thus requires some kind of func- 
tionalization, e.g. through doping or decoration of the CNT 
sidewall [13-21]. Ideally, this type of functionalization could 
be used to control not only the reactivity of the CNT but also 
the selectivity towards specific chemical species. 

In this work we consider the possibility of using CNTs 
doped by 3d transition metal atoms for chemical gas sens- 
ing. We use computational screening to systematically iden- 
tify the most promising dopant candidates for detection of 
three different target molecules (CO, NH3, H2S) under typi- 
cal atmospheric conditions. The screening procedure is based 
on the calculation of two microscopic descriptors: the bind- 
ing energy and scattering resistance of the molecules when 
adsorbed on a doped CNT. These two quantities give a good 
indication of the gas coverage and impact on the resistance. 
For the most promising candidates we then employ a simple 
thermodynamic model of the CNT sensor. In this model, the 
binding energies are used to obtain the fractional coverage of 
the metallic sites as a function of the target molecule concen- 
tration under ambient conditions. Under the assumption of 
transport in the diffusive rather than localization regime, the 



change in CNT resistivity may then be obtained from the cal- 
culated coverages and single impurity conductances. 

We find that oxidation of the active metal site passivates 
the sensor in the case of doping by Ti, V, Cr, and Mn un- 
der standard conditions (room temperature and 1 bar of pres- 
sure). Among the remaining metals, we identify Ni as is the 
most promising candidate for CO detection. For this system 
the change in resistance per active site is generally significant 
(>1 17) for small changes in CO concentration in the relevant 
range of around 0.1-10 ppm. Our approach is quite general 
and is directly applicable to other nanostructures than CNTs, 
other functionalizations than metal doping, and other back- 
grounds than atmospheric air. 

All total energy calculations and structure optimizations 
have been performed with the real-space density functional 
theory (DFT) code GPAW [22] which is based on the projector 
augmented wave method. We use a grid spacing of 0.2 A for 
representing the density and wave functions and the PBE ex- 
change correlation functional [23]. Transport calculations for 
the optimized structures have been performed using the non- 
equilibrium Green's function method [24] with an electronic 
Hamiltonian obtained from the SIESTA code [25] in a dou- 
ble zeta polarized (DZP) basis set. Spin polarization has been 
taken into account in all calculations. 

Metallic doping of a (6,6) CNT has been modeled in a su- 
percell containing six repeated minimal unit cells along the 
CNT axis (dimensions: 15 Ax 15 Ax 14.622 A). For this size 
of supercell a T-point sampling of the Brillouin zone was 
found to be sufficient. The formation energy for creating a 
vacancy (VC) occupied by a transition metal atom (M) was 
calculated using the relation 

£ form [M@VC] = E[M@VC] + nE[C] - £[M@NT] (1) 

where £?[M@VC] is the total energy of a transition metal 
atom occupying a vacancy in the nanotube, n is the number 
of carbon atoms removed to form the vacancy, E[C] is the en- 
ergy per carbon atom in a pristine nanotube, and _E[M@NT] 



Carbon Nanotube Axis 



(a) Adsorption Energy [eV] (b) Conductance Change [G ] 
-2.0 -1.5 -1.0 -0.5 0.0 -1.0 -0.5 0.0 +0.5 +1.0 



IMonovacancy 



□ Divacancy I 



□ Divacancy II 




Empty Monovacancy 
Empty Divacancy I 




V Cr Mn Fe Co Ni Cu Zn 



FIG. 1: Structural schematics and formation energy for a 3d tran- 
sition metal occupied monovacancy (black), divacancy I (gray), or 
divacancy II (white) in a (6,6) carbon nanotube. Formation energies 
of the empty vacancies are indicated by dashed lines. 



is the total energy of the pristine nanotube with a physisorbed 
transition metal atom. We have considered the monovacancy 
and two divacancies shown in Fig. 1. The energy required to 
form an empty vacancy is obtained from 



E lom [VC] = E[VC] + nE[C] - £[NT], 



(2) 



where E[WC] is the total energy of the nanotube with a va- 
cancy of n atoms. 

The calculated formation energies for the 3d transition met- 
als are shown in Fig. 1. From the horizontal lines we see that 
both divacancies are more stable than the monovacancy. This 
may be attributed to the presence of a two-fold coordinated C 
atom in the monovacancy, while all C atoms remain three-fold 
coordinated in the divacancies. When a transition metal atom 
occupies a vacancy, the strongest bonding to the C atoms is 
through its d orbitals [26]. For this reason, Cu and Zn, which 
both have filled d-bands, are rather unstable in the CNT For 
the remaining metals, adsorption in the monovacancies leads 
to quite stable structures. This is because the three-fold coor- 
dination of the C atoms and the CNT's hexagonal structure are 
recovered when the metal atom is inserted. On the other hand, 
metal adsorption in divacancies is slightly less stable because 
of the resulting pentagon defects, see upper panel in Fig. 1 . A 
similar behaviour has been reported by Krasheninnikov et al. 
for transition metal atoms in graphene [21]. 

The adsorption energies for N2, O2, H2O, CO, NH3, and 
H2S on the metallic site of the doped (6,6) CNTs are shown in 
Fig. 2(a). The adsorption energy of a molecule X is defined 
by 

E Rds [X@M@YC] = E[X@M@YC] - E[X] - E[M@\C], 

(3) 
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FIG. 2: Calculated (a) adsorption energy E^s in eV and (b) change in 
conductance AG in units of Go =2e 2 /h for N 2 , O2, H 2 0, CO, NH 3 , 
and H2S on 3d transition metals occupying a monovacancy (top), 
divacancy I (middle), and divacancy II (bottom) in a (6,6) carbon 
nanotube. 



where E[X @M@VC] is the total energy of molecule X on 
a transition metal atom occupying a vacancy, and E[X] is the 
gas phase energy of the molecule. 

From the adsorption energies plotted in Fig. 2(a), we see 
that the earlier transition metals tend to bind the adsorbates 
stronger than the late transition metals. The latest metals in 
the series (Cu and Zn) bind adsorbates rather weakly in the 
divacancy structures. We also note that O2 binds significantly 
stronger than any of the three target molecules on Ti, V, Cr, 
and Mn (except for Cr in divacancy I where H2S is found to 
dissociate). Active sites containing these metals are therefore 
expected to be completely passivated if oxygen is present in 
the background. Further, we find H2O is rather weakly bound 
to most of the active sites. This ensures that these types of 
sensors are robust against changes in humidity. 

In thermodynamic equilibrium [27], the coverage of the ac- 
tive sites follows from 



Q[X] 



K[X]C[X] 



i + E y *W[y] : 



(4) 



where K = k+/k- is the ratio of forward and backward rate 
constants for the adsorption reaction, 



K[X] = exp 



E ads [X}+TS[X] 



(5) 



In these expressions C[X] is the concentration of species X, 
S[X] is its gas phase entropy and T is the temperature. Ex- 
perimental values for the gas phase entropies have been taken 
fromRef. [28]. 
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FIG. 3: Fractional coverage O in thermal equilibrium of Ni in a (a) 
monovacancy, (b) divacancy I, (c) divacancy II and (d) change in 
resistance AR per dopant site as a function of CO concentration in 
a background of air at room temperature and 1 bar of pressure. The 
reference concentration of CO is taken to be Co =0. 1 ppm. Note the 
change from linear to log scale on the y-axis at AR =10 SI. 



For a given background composition we may thus estimate 
the fractional coverages for each available adsorbate for a 
given type of doping. As an example, Fig. 3(a)-(c) shows the 
fractional coverage of a Ni atom occupying a monovacancy, 
divacancy I, and divacancy II, versus CO concentration in a 
background of air at room temperature and 1 bar of pressure. 
Due to the relatively small binding energy of N2 and H2O as 
compared to O2 and CO, all Ni sites will be either empty or 
occupied by O2 or CO. In particular, Ni in a monovacancy 
(top panel of Fig. 3) will be completely oxidized for all rel- 
evant CO concentrations. For the Ni occupied divacancy II 
structures we find the coverage of CO changes significantly 
around toxic concentrations (~10 ppm). 

To estimate the effect of adsorbates on the electrical con- 
ductance of doped CNTs, we first consider the change in con- 
ductance when a single molecule is adsorbed on a metal site of 
an otherwise pristine CNT. In Fig. 2(b) we show the calculated 
change in conductance relative to the metal site with no ad- 
sorbate. In contrast to the binding energies, there are no clear 
trends in the conductances. The sensitivity of the conductance 
is perhaps most clearly demonstrated by the absence of cor- 
relation between different types of vacancies, i.e. between the 
three panels in Fig. 2(b). Close to the Fermi level, the conduc- 
tance of a perfect armchair CNT equals 2Gq. The presence 
of the metal dopant leads to several dips in the transmission 
function known as Fano antiresonances [20]. The position 
and shape of these dips depend on the rf-levels of the transi- 
tion metal atom, the character of its bonding to the CNT, and 
is further affected by the presence of the adsorbate molecule. 
The coupling of all these factors is very complex and makes 
it difficult to estimate or rationalize the value of the conduc- 
tance. For the spin polarized cases, we use the spin-averaged 



conductances, i.e. G = (G^ + G\)/2. 

Next, we estimate the resistance of a CNT containing sev- 
eral impurities (a specific metal dopant with different molecu- 
lar adsorbates). Under the assumption that the electron phase- 
coherence length, I,/,, is smaller than the average distance be- 
tween the dopants, d, we may neglect quantum interference 
and obtain the total resistance by adding the scattering resis- 
tances due to each impurity separately. The scattering resis- 
tance due to a single impurity is given by 



R S (X) = l/G(X) - 1/(2G ) 



(6) 



where G(X) is the Landauer conductance of the pristine CNT 
with a single metal dopant occupied by molecule X and 
1 /(2Gq) is the contact resistance of a (6,6) CNT. 

We may now obtain the total resistance per dopant site rel- 
ative to the reference background signal as a function of the 
target molecule concentration 



AR 



^R s (X)(e[X,C]-Q[X,C ]), 



(7) 



where N is the number of dopants, Q[X, C] is the fractional 
coverage of species X at concentration C of the target and Cq 
is the reference concentration. Notice that the contact resis- 
tance drops out as we evaluate a change in resistance. 

In Fig. 3(d) we show the change in resistance calculated 
from Eq. (7) as a function of CO concentration for Ni occu- 
pying the three types of vacancies. The background reference 
concentration of CO is taken to be Co = 0.1 ppm. For the 
monovacancy there is very little change in resistivity. This is 
because most active sites are blocked by O2 at relevant CO 
concentrations, as shown in the upper panel of Fig. 3. For Ni 
in the divacancies there is, however, a change in resistance on 
the order of IS! per site. For concentrations above ~1 ppm, 
the CO coverage of Ni in the divacancy II increases dramati- 
cally and this leads to a significant increase in resistance. 

We now return to the discussion of the validity of Eq. (7). 
As mentioned, the series coupling of individual scatterers 
should be valid when 1$ < d. However, even for 1$ > d 
and assuming that the Anderson localization length, Zi oc in 
the system exceeds 1$, Eq. (7) remains valid if one replaces 
the actual resistance R by the sample averaged resistance (R) 
[29]. At room temperature under ambient conditions, interac- 
tions with external degrees of freedom such as internal CNT 
phonons and vibrational modes of the adsorbed molecules 
would rapidly randomize the phase of the electrons. There- 
fore Eq. (7) should certainly be valid in the limit of low dop- 
ing concentrations. On the other hand, the total number of 
dopants, N, should be large enough for the statistical treat- 
ment of the coverage to hold. Finally, we stress that Eq. (7) 
represents a conservative estimate of the change in resistance. 
In fact, in the regime where 1$ > k oc , i.e. in the Anderson 
localization regime, the resistance would be highly sensitive 
to changes in the fractional coverage of active sites. Calcula- 
tion of the actual resistance of the CNT in this regime would, 
however, involve a full transport calculation in the presence of 
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all N impurities. At this point it suffices to see that the con- 
servative estimates obtained from Eq. (7) predict measurable 
signals in response to small changes in concentration of the 
target molecules. 

To our knowledge, controlled doping of CNTs with transi- 
tion metal atoms has so far not been achieved. It has, how- 
ever, been found that metal atoms incorporated into the CNT 
lattice during catalytic growth are afterwards very difficult to 
remove [30]. Furthermore, it has been shown that CNT vacan- 
cies, which are needed for the metallic doping, may be formed 
in a controlled way by irradiation by Ar ions [31]. This sug- 
gests that metallic doping of CNTs should be possible. 

In summary, we have presented a general model of nanos- 
tructured chemical sensors which takes the adsorption en- 
ergies of the relevant chemical species and their individual 
scattering resistances as the only input. On the basis of this 
model we have performed a computational screening of tran- 
sition metal doped CNTs, and found that Ni-doped CNTs are 
promising candidates for detecting CO in a background of air. 
The model may be applied straightforwardly to other nanos- 
tructures than CNTs, other functionalizations than metal dop- 
ing and other gas compositions than air. 
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